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ABSTRACT 



The theory of the ionization and radiative corrections to 
electron scattering cross sections has been reviewed to determine 
the formulas most valid for the experimental arrangement at the 
NPGLINAC. Experimentally it has been determined that the cor- 
rections most nearly account for undetected electrons when the 
lower limit of the measured scattering spectrum is at least four 
half-widths below the peak of the spectrum. It has also been 
found that the corrections, most likely the Bethe-Heitler cor- 
rection, induce a positive error of several percent to the cross- 
section when thick targets are used. 
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I. INTRODUCTION 



Linear electron accelerators (LINAC) are well suited for the 
study of nuclear structure. The interaction between a nucleus and 
a high energy incident electron, which is entirely electromagnetic, 
is thought to be well understood. By momentum analyzing electrons 
of initial energy which are scattered through an angle 9 into a 
solid angle d Cl> one can determine the differential cross section 
for scattering, dcr/d^. This cross section can yield much infor- 
mation about the charge structure of the target nuclei. 

Hofstadter [l,2] has written a full review of the theory and tech- 
niques of elastic electron scattering. 

It is well known that these high energy electrons radiate a 
portion of their energy (br emsstrahlung ) in the presence of nuclei 
or other electrons. This br emsstrahlung leads to a broadening of 
the scattered electron distribution such that there is a finite 
probability of scattered electrons with energy approaching zero. 

A broadening of the distribution is also due to energy loss by 
ionization. Since in most experiments it is not feasible to mea- 
sure the spectrum down to zero energy, there is a lower energy 
limit to the measured spectrum. It is therefore necessary to cal- 
culate and apply corrections to the experimentally determined cross 
sections to account for the undetected low-energy scattered 
electr one . 

As the precision and accuracy of scattering experiments 
improve, these corrections become more and moru significant. If, 
for example, a ten percent correction were only known to ten 
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percent, this would result in only a one percent contribution to 
the uncertainty in the experimental cross section. But obviously 
such uncertainties are not acceptable when attempting to do an 
experiment to one percent accuracy. 

This paper is concerned with three distinct corrections. The 
first of these corrects for energy loss due to Landau straggling 
(ionization). The other two are properly called radiative cor- 
rections. The Schwinger correction compensates for br emsstrahlung 
in the field of the target nucleus (scattering center) while the 
Bethe-Heit ler correction deals with br emsstrahlung in the field 
of atomic electrons and other nuclei in the target. The cross 
section for the first process is proportional to the target thick- 
ness while for the lattec it is proportional to the square of the 
thickness . 

The following sections include a review of the theory for each 
of the three processes, as best suited for application at the Naval 
Postgraduate School LINAC (NPGLINAC) , that is, for scattering ex- 
periments in the 30 to 100 MeV range where the scattered electrons 
are detected and momentum analyzed while the recoil nuclei are not 
detected. The NPGLINAC is described in detail in theses by Barnett 
and Cunneen [ 3 ] and Midgarden [Zf]. 

It should be noted that in actual applications the corrections 
are multiplicative corrections to the experimentally determined 
cross sections. However, to be consistent with most publications 
in this field, the sections concerned with the theory use the form 



dfT J 



/exp 



da \ 

■dii ,k --- s -* 



i ”s 
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/d a \ . 

where K^, S g , are the corrections, / 1S the corrected cross 

section while ^ d' H " ) exp '*' s t * le GX P er i men tally determined raw 
cross section. 

Since the Bethe-Heitler correction is proportional to the 
square of the target thickness, its validity can be checked ex- 
perimentally by measuring on targets of different thickness. 

Such experiments were performed and the results are reported. The 
proper applications of the corrections, for example, the choice of 
lower energy limit of the spectrum, have also been examined. 
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IONIZATION CORRECTION 



II. 



As a monoener get ic beam of electrons passes through a scat- 
tering target the resulting spectrum of scattered electrons is 
broadened, in part, by ionization loss, or Landau straggling. 

This straggling induces an error in the cross-section determination 
because the scattered spectrum is not normally measured to zero 
energy. Thus a correction should be applied to account for these 
undetected electrons. This correction has largely been ignored 
at the NPGLINAC and other accelerators. For targets and energies 
encountered at the NPGLINAC, the ionization correction is typi- 
cally about one percent. Until recently the accuracy of cross- 
section measurements was insufficient to warrant the inclusion of 
such small corrections. Since many cross-section measurements are 
made relative to another nucleus, the correction is even less 
important. A review of this correction is provided for complete- 
ness and possible future use. 

Since in electron scattering experiments it is neither desirable 
nor feasible to measure the scattered spectrum down to zero energy, 
there is a lower energy limit of the detected spectrum, E \ , so that 
AE = E 0 - E* . This is illustrated in Fig. 1. Thus a fraction of 
the spectrum, V(AE), is lost by ionization. Then the ionization 
correction can be * xpccted to have the form: 



K i = 1 - V( AE) . 



This is applied to the exper imental cross-section: 
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Landau [ 5 ] first calculated the energy distribution of 

scattered electrons when monoener getic electrons are incident on 

2 

a target of thickness t(g/cm ). Following this theory, the 
fraction lost can be calculated by integration of the Landau 
function 0( AE) , 

V l (AE) = 1 - Y ( AE) 

* E 

o 

0 ( AE)dE. 
o 



Y ( AE) 






0 ( AE) dE 



A graph of V^(AE) is given by Landau and also by Symon [6] and is 
reproduced as Fig. 2. 

V (AE) has been approximated by Isabelle and Bishop [ 7 ] who 
have derived an expression for the correction, K , 



K 



L 



1 



1.21 

AE 




( 2 - 1 ) 



where a = 0 . 1537 • Z*p/A* p 2 MeV cm^g 1 . 

The above equation has regions of applicability which must be 
observed : 

a. When AE is of the order of a*t or smaller, the 
approximation of V^(AE) is unsatisfactory. 

b. For very thin targets the Landau theory is only a 
zero-order approximation. This treatment neglects 
fluctuations due to distant collisions, in which the 
atomic electrons cannot be treated as free electrons [8]. 

To improve the Landau correction, Breuer [ 9 ] used V (AE) as 
given by Landau rather than the Bishop- Isabelle approximation. He 
also incorporated the theory of Blunk and Ledsegang [lO] since this 
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Fig. 1. 

Typical elastic spectrum 





Fig. 2. V as given by Landau. 

L* 
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Counts 



shows better agreement with experiments than the Landau theory. 

Since the calculations of Blunk and Leisegang introduce only a 

broadening of the distribution, the Landau V^(AE) can be used by 

substituting the Blunk and Leisegang values for the full width at 

half maximum ( half-width) , Breuer has plotted this half-width 

2 

as a function of b , where 

b 2 = q Q z Z * /3 /(a-t) 
q = 20 eV 

Q = Q-jPt = mean energy loss for thickness t 

— 2 

- mean energy loss per g/cm . 

2 

The plot of F /(a-t) as a function of b is given in Fig. 3* 
can be calculated by an empirical equation given by 
Sternheimer [ ll] : 

Q x = t [b+ 0.43 E x -p 2 + C -a (X-log p/mc) m ]MeV. (2-2) 

P 

The empirical constants for this equation, as determined by 
Sternheimer, are given in Table 1 for selected targets. 

The abscissa values of Fig. 3 are defined as T) , 

V L = AE/1.5l*a-t . 

The Landau relation for the half-width, 

P L = 4 ' a ' t! 

can be substituted into the above relation to give 
V L = AB/( 0.379 r L ). 
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• J 



16 



Fig. 3* Half-width T as predicted by Blunk and Leisegang. 
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Now the Landau half-width should be replaced by the larger and 

more correct Blunk and Leisegang half-width, T : 

J3 

V n = 2.65-AE/r 

o 

For relativistic energies and r? > 2 , v(AE) falls off inversely 

r5 

with AE. From Fig. 2 it is seen that V(rj = 10) = 0.08, and thus 



1 - v(Tj = 10) = .92 . 



Hence the ionization correction can finally be written as 

0.8 



K. = i - 
x T? 



B 



0.8 

TZ3 AE * 



(2-3) 



To apply the ionization correction it is first necessary to 
2 

determine b for the particular experiment. is then determined 

from Fig. 3 and substituted into the expression for . 

2 

In most experiments at the NPGLINAC, b -4 0.5 so that the 
Landau relation can be used: 
r B /(a-t)= 4 . 

Th en 

t? b = . 6625 * AE/ ( a * t ) , 



so that 



K. 

1 



0.8 (a • t ) 
.6625 AE 



(2-4) 
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III. RADIATIVE CORRECTIONS 



The radiative corrections to elastic scattering are of two 
types. The Schwinger correction, called simply the radiative 
correction by some authors, accounts for the emission and re- 
absorption of virtual photons in the field of the scattering 
nucleus, nuclear br emss trahlung , as well as the emission of soft 
photons of energy less than a specified cutoff energy, * The 
Bethe-Heitler correction, also called the radiation tail, accounts 
for the emission of real photons, by br emss trahlung in the field 
of atomic electrons and in the field of nuclei other than the 
target nucleus. 



A. SCHWINGER CORRECTION 

There are sever 1 version!' of the Schwinger correction 
of which are improvements on the original correction first 
by Schwinger [ 12] . His calculations include the following 
assumpt ions : 

a. stationary nucleus 

b. Coulomb point field 

c. first Born approximation 

d. one-quantum emission 

e . p = v/c ^ 1 
with the restrictions 



a . AE < < E^ - m c 

9 2 

b. E^ sin — / me >>1. 



all 

given 
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The Feynman diagrams corresponding to Schwinger f s calculation are 



given in Fig. 4- This correction, b , is applied as 

O 



dCT 

dTF 



exp 



da 

dir 



1-6. 



wher e 

6 _ 



4# 

TT 



(3n 



AE 



13. ,, 2p o . 9 

T2>< fa T" Sln 2 



1.. 17. 1 . 2 < 

2> + TT 2 S1 " ; 



This can be put into a more convenient form [ 13 ]. 

b s ‘f (< fc f . 

wher e 

^ / r\ \ n r . 2 0 x i, , 29 x . 0 

f(9) - Bm[sx n — )0ft(cos — ) + 5 (sin — ) . 

$ is the Spence function^ 14] 



*(*) 



= j x _ %U 



-y| 



dy 



To simplify calculations this is usually seen as 

b s=f (< to A-T !) • 



$(x) 
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$(x) 
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if |x|^ 



if x > 1 
if x < - 



f(«) 



(3-1) 
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Fig, 4. Feynman diagrams corresponding to the 
original Schwinger correction. 
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Fig. 5. Feynman diagrams corresponding to the 

formulas of Meister and yennie and Tsai. 
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Since this calculation was made for potential scattering, the 
energy of the final and the incident electrons are the same, and the 
formula is ambiguous as to the energies to use and the definition 
of AE. If we neglect all energy losses except recoil, the most 
probable energy of the elastically scattered electrons is 

E 3 = E 1 /T] (3-2) 

where n is the recoil factor 

V = 1 + E 1 M -1 (l-cos ©). ( 3 - 3 ) 



It can be shown that the maximum energy of a photon emitted 

in the direction of the scattered electron is AE, while in the 

2 

direction of the incident electron the maximum energy is T) AE. 
To account for this Tsai [ 1 3 3 gives an improved formula: 




Note that when AE "• 0 all of the above equations are infinite, 

although we would expect the measured cross section to go to zero 

as AE 0. This is apparently due to the fact that multiple 

photon emissions have been neglected. Hence Schwinger suggested 

" b S r , 

that (1 - b ) be replaced by e .It was later proved L 1 5 J that 
the term 

should be exponentiated. It is still uncertain [ 1 6 3 if the other 
terms hould be exponentiated, but the difference is negligible 
at energies obtainable at the NPGLINAC. 



2 a 

TT 






(?) 



v AE' 
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Meister and Yennie [l6] give an improved version of the 



Schwinger cor r ection : 




where p is the step function: p(x) = x)Q( 1-x ) . 

This calculation is described by the Feynman diagrams in 
Fig. 5 • The detailed assumption and restrictions for this 
calculation are given b Breuer [ 17 ], and are similar to those 
of the Schwinger calculation. 

The latest and most accurate version of the Schwinger cor- 
rection is given by Tsai [ 18] . This is essentially the Meister 
and Yennie formula without approximations for the Spence functions. 
This formula takes into account both target recoil and a dynamical 
effect due to photon emission. Two-photon emission still has not 
been sufficiently extracted except in order to cancel infrared 
divergence in real photon emission. The assumptions, restrictions 
and Feynman diagrams are the same for the Tsai formulation as for 
the Meister and Yennie formula. Tsai’s expression for the 
correction is: 
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B. BETHE-HEITLER CORRECTION 

The Bethe-Heit ler correction accounts for energy loss by 
br emsstr ahl ung in the field of nuclei other than the target 
nucleus, as well as br emsstr ahlung in the field of atomic 
electrons. A review of br emsstrahlung , screening, and the 
resulting cross-section formulas is given in Appendix A. 

Bethe and Heitler [l9,20,2l] have shown that if the cross- 
section for br emsstr ahlung were given by 

TU =toFx— /E)]- 1 , (3-7) 

o 
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the probability an electron of initial energy E^ being in an 
energy interval dE at E, after traveling thrpugh a target of 



thickness T (in units of radiation length), would be 

E^CNE./B )] 1 '®" 2 

I (E E T) = — 

e 1 F ( T/&n 2 ) 



(3-8) 



Integration of this formula, with the assumption that T is small, 
leads to the original Bethe-Heit ler correction, first given by 
Hofstadter as 



b_ = 



B &7l2 



On 



E 

AE 



To include the effect of target recoil, this should be written 



b 



B 



Qm. 2 



Qm, 




(3-9) 



Betho and Heitler admit that Eq. (3-7) is a poor approximation 
to the correct cross section and was chosen for its simplicity. 

To improve uj )n the correction it is necessary to use a more nearly 
correct expression for the cross-section. However, it is also 
necessary to be able to perform the integrations. 

Tsai [l8] has developed a unique method for the case of 
complete screening. By combining Eq. (3-7) and Eq. (3-8) for 
the case of small T, Eq. (3-8) becomes 



I (E. , E,T) = /- X 

e' 1 ' \ A o 



I ' 



dcr 

dE > 



Q/n 



E x \ T/to2 



(3-10) 



T/h 2 

Tsai shows tha < the term (2n E^/E) ” is a correction for 

multiple scattering which is insensitive to the choice of the 
cross-section formula. Hence it would seem appropriate to replace 
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dCT/dEin Eq. (3-10) by a correct expression, thereby immediately 
obtaining the intensity I . By adding the exact cross sections 
for nuclear and atomic electron br emss trahlung , the total cross 
section for one-photon emission and complete screening can be 
written (with a cross term of < 0 . 5% deleted) 



da 

dE 




elec 



.-1 A 4 

o N 3 





1 + 



E Z+l 
9E 1 Z+§ 



^(183Z 




y 



where X q is the radiation length given by Bethe and Ashkin [ 22 ] 



1_ 

X 




ar Q Z(z + 5 ) 0n(l83z 



Thus Eq. (3-8) becomes 



I e (E 1 ,E,T)= bT(E 1 -E) 




( 3 - 12 ) 



where to a very 




good approximation 



[(Z+ l)/(Z+§) 



M183Z 





(3-13) 



and § is the Wheeler “Lamb factor. 

From this, Tsai evaluates the correction under the assumption 

that the scattering event occurs midway through the target. He 

also incorporates energy loss in an entrance window of thickness 

T. and a final window of thickness 1' to derive the formula 
xw fw 






■{ 



b T. + X bT 
w iw 2 



&n(E /ri LE) 



b T + 
w fw 



\ bT 



j &?(E 3 /AE) J. 



(3-i4) 
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In examining, for example, the proton scattering in a solid target 
such as (CH^)^ br emss tr ahlung from the carbon in the target can be 
accounted for by substituting the effective thickness of the carbon 
for the window thickness. 

The expression for b as derived by Tsai is slightly more 

B 

accurate than Eq . (3-9) and certainly more versatile. It is, 

however, derived on the assumption of complete screening, an 

approximation which begins to fail at energies attainable by 

the NPGLINAC. Unfortunately, the magnitude of the errors induced 

is, as yet, unknown. A more correct theory would require use of 

a cross-section formula which includes the effects of "inter mediate" 

screening. These formulas, which include atomic form factors, are 

not known explicitly and determination of the correction would 

require a double numerical integration, with no hope of an ex- 

2 

plicit expression for the Bethe-Heitler correction . There was 
insufficient time to pursue this idea further for inclusion in 
this paper . 



'Incorporating the effects due to screening also means 
re-evaluating the Wheeler -Lamb factor which in turn implies 
changes in the radial on length. 
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IV. CHOICE OF AE 



The shape of an elastic spectrum near the peak is largely due 
to the energy spread of the incident electron beam. A broadening 
of the peak is also due to the finite width of the entrance slit 
as well as the ionization and radiative corrections. The cor- 
rections should be able to account for electrons lost by ion- 
ization and by radiation for any reasonable choice of AE. However 
AE must be chosen sufficiently large to include the energy spread 
due to the two machine effects. 

Tsai [ 13 ] gives requirements on the choice of AE. If the 
incident beam has an energy spread of AE^ , the equation 

E 3 = Ej/17 

is used to calculate the energy spread of the scattered electrons: 




Tl is the recoil factor given in Eq. (3“3). Similarly Tsai shows 
that the energy spread due to the finite width of the entrance 
slit is 




Hence AE should be chosen such that 

-2 

AE>>AErf , 

and 

2 

AE > (E 3 / M) sin 9 a9. 
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In the development of the Tsai formulation of the Schwinger 
correction there is an assumption 

E( 1 + 2E 1 /M)<< E . 

This puts n upper limit on the choice of AE, but is not a 
severe restriction. 
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V. EXPERIMENTAL CONSIDERATIONS 



A. COUNTING SYSTEM 



The design and operation of the NPGLINAC has been described 
in detail in several theses [3>43. The most significant recent 
modification has been the installation of a ten-channel counting 
system. Ten plastic scintillators have been positioned vertically 
in the focal plane of the l6 inch double-focusing spectrometer. 

The spectrometer is described in a thesis by Oberdier [ 23 ]. 

Electrons which enter the spectrometer with an energy (more 
properly, momentum) corresponding to a given spectrometer setting 
follow a path through the spectrometer such that they focus at 
the intersection of the focal plane and an extension of the 
central path of the spectrometer. Electrons of greater energy 
are focused above this point in the focal plane while less 
energetic electrons are focused at a lower point. The distance 
from the central point is given roughly by the dispersion formula 



wher e 




(5-1) 



r = central radius of the spectrometer = 16 inches 
D = dispersion constant = 3*92 

= momentum difference from central focusing point. 




The ten scintillators are arranged above and below the 
central point (channel five very nearly corresponds to the 
central energy) so that each one "sees" a different energy. 
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Immediately behind the front counters is a large vertical backing 
scintillator. To minimize the detected background due to the 
intense radiation field in the accelerator end station, only 
electrons which pass through both one of the front counters and 
the back counter in coincidence are registered as scattered 
electrons. The coincidence system is illustrated schematically 
in Fig. 6. 

To get a reasonable density of points defining a scattering 
peak it is necessary to change the central energy of the spec- 
trometer in small steps of one-third to one-sixth the average 
counter separation (0.10 to 0.05 MeV at = 9° MeV), as well as 
one or more large steps ( ~2.5 MeV at E 1 = 90 MeV) to define the 
radiation tail. This requires a rather involved data unfolding 
procedure before a usable spectrum is produced. We follow 
closely the procedure given by Suclzle [ 24 , 25 ]. The energy of 
each « ounter as a function of the central energy has been deter- 
mined by plotting, at several incident energies, the elastic 
peak as it occurs in each channel as a function of the spectro- 
meter setting. The resolution of the i-th channel is defined as 



E. 



1 




E 



i-1 



)• 



The efficiency of each channel is determined for each run by 
taking at least five overlapping data sets on a smooth portion of 
the radiation tail and fitting these points to a second order poly- 
nomial. The efficiency o i each channel is then found from the 
difference of .he data points and the polynomial fit. 
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Pig. 6. Counting System 
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B. EXPERIMENTAL METHOD 



1 . Basis of the Experiment 

The purpose of this experiment is to examine the adequacy 

of the radiative corrections as a function of target thickness, and 

the importance of the cutoff energy, E , to the corrections. This 

3 

leads to the straight-forward procedure of determining the cross 
section for several thicknesses of the same target material at a 
given incident energy and scattering angle. Since we are only 
interested in the relative differences between targets, it is 
sufficient to determine the area under the peaks. If the cor- 
rections are correct, the ratio of the corrected area to the 
target thickness will be constant. Similarly, for any single peak, 
the corrected area should be independent of AE, provided the 
conditions discussed in Section IV are met. 

These corrections have been examined for six thicknesses 
of graphite, p = 1.44l g/ cm (.028 to .250 inches) and seven 
thicknesses of aluminum, p = 2.543 g/cnr (.0104 to .1354 inches). 
Incident energies of 35, 55 and 90 MeV were used. The scattering 
angle was 90 degrees in all runs. The dependence of area on AE 
was examined on several targets at different energies. 

2 . Internal Consiste n cy of the Data 

In this work it was necessary to measure the area 
(counts) of an elastic peak but not necessary to actually de- 
termine the cross section. Therefore errors in the measurement 
of the absolute valu.es of experimental constants (incident energy, 
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scattering angle, target angle, solid angle, integrator capacitance, 
SEM efficiency and others) do not enter into these relative experi- 
ments. However, it is important to insure that none of these fac- 
tors vary during a data run. 

To avoid effects which may be due to the tuning of the 
accelerator, such as variations in incident energy, energy reso- 
lution of the beam, and spot size of the beam, it has not been 
attempted to compare different runs. A run, consisting of a 
spectrum from each target, was completed without retuning the beam. 
It is felt that the deflection magnets, which define the incident 
energy, drift less than 0.1% over several hours, while short term 
variations are undetectable. 

To avoid errors due to hysteresis in the sp .ctrometer , the 
magnet current, and hence the central energy, is monotonica lly 
decreased during the run to span the peaks. Because of drift in 
the current regulated supply or temperature change in the magnet, 
the field of the magnet is not exactly proportional to the current 
in the windings. Hence the field is sensed by a r otat ing-coil 
fluxmeter , and maintained at a specified value by a balancing 
circuit. The field is held constant to approximately ^.005 MeV 
in the central energy. The precision potentiometer used in the 
balancing circuit to define the field value is linear to 0.1% 
throughout its range. 

It is critical that the SEM efficiency and its associated 
beam current integrator be constant over the time of the experi- 
ment. Since the beam intensity must vary during a run to account 
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for the difference in target thickness, the efficiency of the 
SEM must be independent of beam current, or a correction must be 
known. This is fully discussed in Appendix B. 

There is some question as to the validity of the efficiency 
determination for the individual counters. To minimize such errors, 
all peaks have been taken so that the top of the peak is defined 
entirely by channel seven. Thus the efficiencies, even if in- 
correct by several percent, are applied similarly to all peaks. 

With the previous one-channel system it was determined 
that a counting rate correction of the form 

C T = G [l .0 + 0.003c/t] 

wher e 

= corrected counts 
C = observed counts 
t = integration time in seconds 

should be applied to the observed counting data. This same cor- 
rection has been applied to data taken on the ten-channel system, 
although it is now apparent that this is too large. However, 
roughly the same counting rate, ten counts per second, was main- 
tained when running each peak, so that errors will again cancel 
out. It is estimated that accidental coincidences do not con- 
tribute more than a 0.1% error to the cross sections. 

The most persistent problem encountered during this experi- 
ment was the detection of an excessive number of electrons in the 
upper channels when counting on a radiation tail. This problem 
was partially solved by simply deleting obviously bad points before 
the numerical integration of the peak. A discussion of this problem 
ill be found in Appendix C. 
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VI. DATA ANALYSIS 



A. DETERMINATION OF CROSS SECTION 

The data from an electron scattering experiment is in the form 
of a counts versus MeV (properly MeV/c) spectrum of the scattered 
electrons. This spectrum is unfolded from the raw data by the 
procedure mentioned in Section V. To find the actual number of 
scattered electrons it is necessary to first divide out the 
momentum acceptance of each counter, Ap . To a very gocJ approxi- 
mation this is given by the dispersion formula 

a/r 



Ap = 



D-p 



(6-1) 



where now a is the effective height of each counter, approxi- 
mately 7/32 inch. The resulting spectrum is thus one of counts 
per MeV versus MeV, so that if the area of the peak is measured 
it has units of counts, the number of scattered electrons. The 
number of scattered electrons per millivolt integration of the 
incident beam is proportional to the cross section. 

Determination of the area was by numerical integration of the 
unfolded spectrum. A three-point method, similar to Simpson T s 
Rule, was chosen. Since the data points were not equally spaced, 
it was necessary to exactly fit each set of three points to a 
parabola, and to explicity perform the integration under the 
curve. It is felt that this method is superior to the commonly 
used trapezoidal rule. 

Once the areas have been determined, it is a simple procedure 
to apply the corrections and normalize the areas to a given tar- 
get thickness. All the corrections except the most complicated 
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of the Schwinger corrections, Eq. (3-6), lend themselves to 
computation by hand or with a calculator . 

The data is presented in two forms. To examine the validity 
of the Bethe-Heitler correction the ratio of the corrected area to 
the target thickness versus target thickness has been plotted for 
each run. Straight line least squares fits have been included. 

To examine the choice of AE, the corrected area versus AE has been 
graphed for several targets at different energies. 

B. EXPERIMENTAL ERROR 

The statistical error in the cross section due to the random 
nature of scattering reactions is very nearly 




where O' is the standard deviation and Sn . is the total number of 

x 

counts in a peak. Typical statistical uncertainties vary from 
one-half of one percent to one percent for the experiments re- 
ported in this paper. 

A significant contribution to the uncertainty of the measured 
cross section is introduced by the numerical integration. By 
adding or deleting points on a peak and repeating the integration, 
this uncertainty has been estimated to be about one percent 
(slightly greater for the thin targets which have fewer points 
defining the peak and slightly less for the thick targets). 
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The target thickness is known to one percent (with smaller 
errors for the thicker targets). For the experiments which com- 
pare the cross section as a function of target thickness, these 
errors propagate to a total uncertainty of about 1.5%. Experiments 
to determine the proper choice of AE are not subject to these 
uncertainties since only one peak is integrated, and the successive 
integrations over the peak are performed only with variations in 
the lower cutoff energy. The relative uncertainty in thes' suc- 
cessive integrations is estimated to be less than one-half of one 
percent . 

Because of the relative nature of all the exper iments , 
systematic errors are thought to be insignificant. 
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VII. RESULTS 



A. CHOICE OF AE 

The data from the experiments to determine the appropriate 
choice of AE are shown in Figs. 7 to 13, and summarized in Table II. 
Figure 7 is a plot of the cross section versus AE for carbon data 
taken on the Mark III linear accelerator at Stanford University [ 26 ]. 
The cross section does not approach an symptotic value even for a 
AE of 20 half-widths. Figures 8 to 13 show typical experimental 
results obtained with the NPGLINAC. In almost all cases an 
syn, totic value was reached for a AE of four half-widths (six 
half-widths in the worst case, E^ = 35-0 MeV) . No significant 
difference was noted when the energy defining slits were opened 
from the normal position of .0920 inches to .1840 inches, al- 
though the resolution ^ 1 of the incident beam, AE^/E^ , changed from 
0.4% to 0.8%. Note that these values are less than two half-widths 
of the elastic peak, so that this is the expected result. 



3 

Cross sections were calculated by the procedure outlined for 
the NPGLINAC; however, Bumiller and Dally [ 26 ] report similar 
results using the same data. 

^AE n /E-. = x/r 
11 o 

x = opening of the slits in inches 

r^^ effective radius of the deflection magnet, 23 inches 
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B. ACCURACY OF THE CORRECTIONS 



The data from the exper iments to determine the validity of 
the corrections to the cross section are shown in plots of cross 
section versus target thickness in Figs. 1 4 to 22, and summarized 
in Table III. In all cases the experimental cross section is con- 
siderably larger for the thicker targets. No systematic errors 
are known which could account for this fact. In fact, the possible 
errors, change in SEM efficiency with beam current (APPENDIX B) or 
multiple scattering events, would tend to magnify the observed 
effect . 
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Summary of Data, Cross Section as a Function of Thickness 
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of half-widths for integration. 
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Fig. 17. 
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VIII. CONCLUSIONS 



The ionization and radiative corrections are most accurately 
given by Equations (2-3), (3-6), (3-14). Thus the cross section as 

determined by a scattering experiment is 



By integrating the experimental spectra for several successively 
lower cut-off energies and applying the corrections, it was deter- 
mined, within experimental error, that a AE of four half-widths is 
more than adequate, even for thick targets (0.01 radiation lengths) 
even with the energy defining slits opened more than usual 
(slits = 200. = .184 inches). The consistency of this data in- 
dicates that the unfolding procedure and the numerical integration 
was done properly. Data taken on the Mark III accelerator does 
not converge even when the integration is carried out to 20 half- 
widths. This discrepancy has not been explained. 

For thick targets these corrections lead to cross sections 
which are too large by several percent. This may be due to the 
fact that the assumption of complete screening in the Bethe- 
Heitler theory begins to break down at energies in use at the 
NPGLINAC . This error has not been estimated. 
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APPENDIX A. BREMSSTRAHLUNG 



A complete theory of br emss tr ahlung has been developed by 
Bethe and Heitler [l9,20,2l]. Other authors have developed 
similar theories, but they lead to essentially the same results. 
Since knowledge of the br emss tr ahlung cross sections is basic to 
the Bethe-Heitler radiation correction, it is appropriate to 
review the basic features of the Bethe-Heitler theory. 

When an electron of initial energy passes through the 
electromagnetic field of a nucleus, the interaction of the two 
fields results in an acceleration to the electron so that it must 
emit radiation. The emission of a quantum of momentum k/c leads 
to a final electron state, of energy E: 

k = hV=E 1 - E . 

The problem is thus to find the transition probability, or cross 
section to the final state E (momentum p) . 

The analysis involves the interaction of the electron with the 
radiation field^, H, and the nuclear field^, V, both of which are 
treated as perturbations. There are thus two intermediate states, 
one caused by H where the quantum k is emitted, so that the electron 
has raomentuu 




and a competing intermediate state caused by V, where the electron 

has momentum ■/ 

~*// 

p = p + K 



H = - e CT A, where A is the vector potential, OL the velocity. 

2 

Assume a Coulomb field, V = Ze /r . 
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and k is emitted with conservation of momentum in the transition 

from intermediate to final states. By appropriately summing over 

the intermediate states one can determine the matrix elements for 

7 

the total transition, . This leads directly to the differential 

g 

cross section . This has been integrated over all angles by 
Bethe [2l] to derive the cross section (differential in energy) 
for the case of no screening: 



da 

dE 



2 2 

Z ar 

o 

[Bj- e] 



4 



1 + 





2E ][ E 

nc 2 [E ]L -E] 




(A-l) 



This result was derived on the assumption of a pure Coulomb 
field, but in certain cases the screening of the nuclear field 
by the atomic electrons necessitates significant changes. The 
charge distribution changes the original potential 

y _ J exp [ i(q-r )/ficJ _ hjr£_c ^2 

q 

, 2 2 

V = [ z . K(q) f 

q 

where F is the atomic form factor of the Fermi-Thomas model 
F(q) =Jp(r)exp (q-r )]dT , 

with p( r) the density of atomic electrons. Hence screening may be 

2 

accounted for by replacing Z in the original differential cross- 
section formula by [z-F(q)l . 



Equation (25-8) in Reference 20. 
Equation (25-13) in Reference 20. 
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For screening to be significant, the field must be screened 
appreciably for "impact parameters" which contribute most to the 



br emss tr ahlung process. F(q) becomes comparable to Z for 

2 

137 



q £ as* z l/3 



which is tic times the reciprocal of the atomic radius. It can 
also be shown that there is a minimum value of q, 

2 2 

- (roc ) h V 

Smin 2E^E 



Then screening is effective if 

1 2 - i 

E 1 E/hV>- 137 me Z 3 

Numerical integration over the angles involved by Bethe leads 
to the general formula 



da 

dE 



2 

z a 



E (E -E) 
o v o ' 





0 2 (y)-| fa z) 



(A-2) 



wher e 



y = 



100 me h V 



E E Z 
o 



1 

3 



0±, 0^ are given in Fig. 23. Y determines the effect of screening. 
If y = 0 the screening is said to be complete while for y > > 1 
screening may be neglected. 

In the case of complete screening, Eq. (A-2) reduces to 



da 

dE 



z 2 < 

Le^eJ 




2 _E_ 

3 E, 



bn 183 Z 




(A-3) 
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Fig. 23. Influence on Screening by Atomic Electrons 
(high energy). & 1 > 0 2 used in Eq. (A-2) . 




(lower energy). C(y) used in Eq. (A-5). 



55 



Two different expressions may be written for the intermediate 



cases : 

y < 2: 



s-S=i‘l(-r; 



1 



2 _E_ 

3 E, 



0 2 (y) ! 



-a 



(a-4) 



J 



2<y<i 5: 



2 2 

Z ar 

d cr o 



dE Le^eJ 



4 







p 


, / E \ 


2 2 E 




I 1 (\ 


l‘35T 





2E 1 E 

2 "* 2 
me hV 



- 4 - c(y) 



(A-5) 



c(y) is plotted in Fig. 24. 

Wheeler and Lamb [ 27 ] have developed an expression for the 
br emsstrahlung cross section for all the electrons in the atom, 
In the limit of complete screening this is 



/'derj 

^ /electr ons 



ZQ!r 4 
o / 

Ib^eJ 4 



1+ 



\ 2 2 _E_ 

) ” 3 E 



Bn (1440 z 3 J + | I • (A-6) 



This leads to the rheeler-Lamb factor 5: 

zd?) 4 - f' 

\ / electr ons ~ £72 y 1440 Z 

/do ) / 1 



I = 



\ /nucleus £ 72(183 Z 3 

a term which occurs in many equations, 
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APPENDIX B. SECONDARY EMISSION MONITORS 



The secondary emission monitor (SEM) was first described by 
Tautfest and Fletcher [28]. An SEM consists of a set of thin 
metal foils placed in the beam so that the electrons are normally 
incident on the foils. Alternate foils are negatively biased to 
several hundred volts, so that secondary electrons ejected from 
the foils by the beam are collected on the positive plates. The 
resultant current, which is proportional to the beam current, can 
be monitored directly with an ammeter or collected on a capacitor 
used in conjunction with vibrating reed electrometer to integrate 
the current. To determi e the efficiency of an SEM, it must be 
compared to an absolute monitor such as a. Faraday cup. Because 
of the high gamma and neutron radiation produced by the Faraday 
cup it cannot be used to monitor the beam during an experiment. 

In most experiments performed at the NPGLINAC, a small, 
three-foil SEM was positioned inside the target chamber immedi- 
ately downstream of the target, with a larger SEM outside the 
target chamber, but under continuous vacuum. The latter is used 
to integrate the beam current, while the small SEM monitored the 
beam current with a Beckman ammeter . 

It was found that this arrangement is not satisfactory for 
experiments involving targets of different thickness. As the 
thickness of the target was increased, the current of each SEM 
was observed to increase (an increase of about 15% was noted for 
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o 

a .2 5 inch carbon target for an incident energy of 35 MeV) . To 
avoid this problem, the small SEM was moved upstream of the 
scattering target and used to integrate the current, while the 
large SEM was used with a Beckman ammeter. 

The upstream SEM, composed of three aluminum foils .0025 
inches thick, is actually a scatter er and served to increase the 
beam spot size. But even with the upstream SEM, and at low 
energies (35 MeV), the spot size never exceeded 3/16 inch in 
diameter. The primary disadvantage of this arrangement was that 
the target chamber , and thus the SEM, was occasionally exposed 

to air. However, the vacuum was attained at least eight hours 

-5 

before beginning the experiment and maintained below 2 x 10 mm Hg . 

It has been reported by several authors [ 29 , 30] that the 
efficiency of these monitors changes significantly when first 
exposed to the beam. A graph of the SEM efficiency as a function 

of time is shown in Fig. 25* On these runs the beam current was 

- 8 

held constant at 8 x 10 amp. The efficiency settled down to a 
constant value 1.5%) after, at most, two hours exposure to the 
beam. 

More important is the SEM efficiency as a function of beam 
current, since the beam current had to be changed considerably 
from thi,, to thick targets during each run. Bumiller and Dally [30] 



This increase is thought to be due to additional secondary 
electrons produced in the target which are collected by the 
SEM. 
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Efficiency > Efficiency 
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have reported a decreasing efficiency of 10% for one of their 
SEMs as the beam current decreased two orders of magnitude. We 
apparently see a reverse dependence on beam current, as shown in 
Fig. 26. 

It is still not clear if a correction should be applied to 
the data to include this effect. This was wholly unexpected, 
and if correct, would tend to increase the slope of the cross 
section versus thickness curves described previously. 
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APPENDIX C. SCATTERING WITHIN THE SPECTROMETER 



The most persistent problem encountered during this experiment 
was a "ghost" of the elastic peak seen on the radiation tail of 
the peaks. The manifestation of this effect is the detection of 
extraneous electrons, primarily in channels eight, nine and ten, 
when the spectrometer is driven at least two percent below the 
energy of the elastic peak. A typical example of this effect is 
the thick target spectrum illustrated in Fig. 27. This problem 
can be minimized by the use of thin targets, where, in most cases, 
the peak can be measured to four half-widths without a major shift 
of the spectrometer . 

It was first thought that this effect -was due to scattering 
through the upper portion of an aluminum flange holding a vacuum 
window to the "snout 11 of the spectrometer. Initially a 1-1/8 inch 
block of lead was placed at the top of the snout to attempt to 
absorb the undesired electrons. This served only to increase the 
overall background. Then the flange was milled to one-half its 
original thickness with sharp edges filed so as not to interfere 
with the beam. This also showed no noticeable effect. The snout 
is illustrated in Fig. 28. 

The rotating coil of the fluxmeter extends about two inches 
into the vacuum chamber of the spectrometer. It was removed for 
a single run to examine its contribution to the "ghosts". No 
change was m ted. 
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Fig. 28 . 

Spectrometer Snout Assembly 
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It is now believed that these extraneous electrons are due 



to a small angle scattering of electrons which comprise the 
elastic peak, as they graze the walls of the snout and a small 
overhang at the top of the spectrometer. Order of magnitude 
calculations can be made of the position of the elastic peak by 
the dispersion formula 

f - =£• • ( 5 - 1 ) 

For 35 MeV electrons to be shifted three inches up the focal 
plane (where they first hit the flange) a shift in the central 
energy of the spectrometer amounting to 



= 4 . 7 %= 1.5 MeV. 

p 3.92 



is required. This corresponds to what is observed. Since this 
effect is observed as the spectrometer is stepped down as much 
as six percent more, it is unlikely that the ghosts are due to 
scattering from a single point, but rather scattering over a 
large region such as the entire top wall of the snout. 

It would seem reasonable to construct a new snout of the 
greatest possible vertical dimension. This would serve to move 
scattering areas as far from the counters as possible. Any 
scattering from the walls of the spectrometer cannot be corrected 
however, lowering the counters would help remove them from the 
worst region. 
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